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Abstract 

Smoldering  combustion  is  the  slow,  low-temperature,  flameless  burning  of  porous  fuels  and  the  most 
persistent  type  of  combustion.  It  is  the  driving  phenomenon  of  wildfires  in  peatlands,  like  those  causing 
haze  episodes  in  Southeast  Asia  and  Northeast  Europe,  but  is  poorly  understood.  In  this  work,  we  develop 
a  comprehensive  1-D  model  of  a  reactive  porous  media,  using  the  open-source  code  Gpyro,  to  investigate 
smoldering  combustion  of  natural  fuels  with  an  emphasis  on  the  roles  of  the  moisture  and  inert  contents. 
The  model  solves  the  species,  momentum,  and  energy  conservation  equations  and  includes  heterogeneous 
chemical  reactions.  A  previously  developed  5-step  reaction  scheme  for  peat,  including  evaporation  of 
water,  is  adopted  to  describe  the  drying,  thermal  and  oxidative  degradation  during  the  smoldering  combus¬ 
tion.  The  model  predicts  the  transient  temperature,  species,  and  reaction  profiles  during  ignition,  spread, 
and  extinction.  The  predicted  smoldering  thresholds  related  to  the  critical  moisture  and  inorganic  contents 
for  ignition  show  a  good  agreement  with  the  experimental  results  in  the  literature  for  a  wide  range  of  peat 
types  and  organic  soils.  The  influences  of  the  kinetic  parameters,  physical  properties,  and  ignition  protocol 
are  investigated.  This  is  the  first  time  that  a  physics-based  model  of  smoldering  peat  fires  is  developed,  thus 
helping  to  understand  this  important  natural  and  widespread  phenomenon. 

©  2014  The  Authors.  Published  by  Elsevier  Inc.  on  behalf  of  The  Combustion  Institute.  This  is  an  open 
access  article  under  the  CC  BY  license  (http://creativecommons.Org/licenses/by/3.0/). 
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1.  Introduction 

Smoldering  combustion  is  the  slow,  low-tem- 
perature,  flameless  burning  of  porous  fuels  and 
the  most  persistent  type  of  combustion  [1],  Smol¬ 
dering  is  the  dominant  phenomenon  in  megafires 
in  natural  deposits  of  peat  which  are  the  largest 
and  longest  burning  fires  on  Earth.  These  fires 


*  Corresponding  author.  Fax:  +44  (0)20  7594  7036. 
E-mail  address:  G.Rein@imperial.ac.uk  (G.  Rein). 


contribute  considerably  to  global  greenhouse  gas 
emissions,  and  result  in  widespread  destruction 
of  ecosystems  and  regional  haze  events  (e.g.  recent 
megafires  in  Southeast  Asia,  North  America  and 
Northeast  Europe)  [2].  It  is  an  emerging  research 
topic  in  climate  change  mitigation  but  is  poorly 
understood.  For  example,  during  the  1997 
extreme  haze  event  in  Southeast  Asia,  peat  fires 
emitted  the  equivalent  to  13^10%  of  the  global 
man-made  greenhouse  gas  emissions  of  that  year 
[3],  Rein  [2]  has  pointed  out  that  the  atmospheric 
release  of  ancient  carbon  from  the  soil  and  the 
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Nomenclature 

A 

pre-exponential  factor 

Greeks 

cp 

heat  capacity 

y 

radiative  conductivity  coefficient 

dp 

characteristic  pore  size 

S 

emissivity 

c 

heat  of  complete  oxidation 

V 

viscosity/stoichiometric  coefficient 

E 

activation  energy 

p 

bulk  density  (mass  concentration) 

h 

enthalpy 

Ps 

solid  density,  p/ps  =  1  —  i// 

hc 

convective  coefficient 

G 

Stefan-Boltzmann  constant 

h 

nm 

mass-transfer  coefficient 

i 

fraction  factor 

k 

thermal  conductivity 

porosity 

K 

permeability 

6j 

reaction  rate 

m" 

mass  flux 

n 

heterogeneous  reaction  order 

Subscripts 

p 

pressure 

0 

initial 

R/R, 

universal/specific  gas  constant 

oc/ oto 

a-char/a-char  oxidation 

S 

particle  surface  area 

l)/l)o 

/?-char//?-char  oxidation 

t 

time 

a 

ash 

T 

temperature 

d/f 

destruction/formation 

Y 

mass  fraction 

dr 

drying 

Z 

distance 

g 

gas 

IC/MC 

inorganic/moisture  content 

i 

condensed  species  number 

j 

gaseous  species  number 

k 

reaction  number 

p/po/pp 

peat/peat  oxidation/peat  pyrolysis 

w 

water 

sensitivity  of  peat  ignition  to  higher  temperatures 
and  drier  soils  could  create  a  positive  feedback 
mechanism  for  climate  change. 

Peat  can  hold  a  wide  range  of  moisture  con¬ 
tents  (MC1),  ranging  from  10%,  under  drought 
conditions,  to  in  excess  of  300%,  under  flooded 
conditions  [4].  Water  represents  a  significant 
energy  sink,  and  furthermore  natural  or  anthro¬ 
pogenic-induced  droughts  are  found  to  be  the 
leading  cause  of  smoldering  megafires  [5,2^ .  There¬ 
fore,  soil  moisture  is  the  single  most  important 
property  governing  the  ignition  and  spread  of 
smoldering  wildfires  [5-7].  The  critical  moisture 
content  (MCt.)  for  initiating  smoldering  of  various 
boreal  peat  has  been  measured  in  the  range  40- 
150%  in  dry  basis  [7,8],  Drier  than  this  threshold, 
peat  becomes  susceptible  to  smoldering.  The  sec¬ 
ond  most  important  property  is  the  soil  inorganic 
content  (IC 1 ).  As  experimentally  found  by  Frand- 
sen  [6,7],  there  is  a  decreasing  quasi-linear  rela¬ 
tionship  between  MCc  and  ICc:  soil  with  a  high 
IC  can  only  be  ignited  at  low  MC.  Mineral  matter 
acts  as  a  heat  sink  but  also  enhances  the  heat 


1  Moisture  content  (MC)  is  defined  in  dry  basis  as  the 
mass  of  water  divided  by  the  mass  of  a  dried  soil  sample, 
expressed  as  %.  Inorganic  content  (IC  <  100%)  is  defined 
in  dry  basis  as  the  mass  of  soil  inorganic  matter 
(minerals)  divided  by  the  mass  of  a  dried  soil  sample, 
expressed  as  %. 


transfer  via  its  higher  heat  conductivity.  After 
moisture  and  inorganic  contents,  other  important 
properties  are  bulk  density,  porosity,  flow  perme¬ 
ability  and  organic  composition  [2], 

The  spread  of  smoldering  fires  is  controlled  by 
heat  and  mass  transfer  processes  in  a  reactive 
porous  media.  The  computational  studies  on 
smoldering  combustion  in  the  literature  have  only 
included  three  fuels:  cellulose  [1,9],  polyurethane 
foam  [1,8]  and  char  [10],  simulated  with  chemical 
schemes  of  different  complexity,  including  1  [10],  3 
[1]  or  5  [8]  steps.  Ohlemiller  [1]  reviewed  the  early 
attempts  on  simulating  smoldering  combustion 
and  provided  the  governing  equations  in  general 
form.  Rein  et  al.  [8]  numerically  solved  the  1-D 
smoldering  combustion  of  polyurethane  foam 
under  forced  flow  with  a  5-step  kinetics,  and  the 
results  were  compared  to  microgravity  experi¬ 
ments  in  both  opposed  and  forward  propagation 
modes.  He  et  al.  [10]  developed  a  1-D  model  to 
solve  the  in-depth  spread  of  smoldering  for  char 
with  1-step  chemistry  and  compared  it  to  the 
experiments.  Previous  studies  have  not  considered 
simulations  of  peat  fires  or  the  drying  process. 

In  this  work,  we  use  a  comprehensive  1-D 
model  based  on  Gpyro  [11]  and  a  previously 
developed  5-step  kinetics  (including  drying)  for 
peat  kinetics  [12]  to  investigate  the  ignition  and 
spread  of  smoldering  in  a  bed  of  peat  and  other 
organic  soils.  Prediction  of  the  smoldering 
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thresholds  related  to  MC„  and  IC£  are  the  empha¬ 
sis.  The  computational  results  are  compared  with 
experiments  [6,7],  The  influences  of  the  kinetic 
parameters,  ignition  protocol,  and  physical  prop¬ 
erties  are  investigated. 


2.  Computational  model 

Frandsen  [6,7]  conducted  two  sets  of  pioneer 
experiments  to  determine  the  smoldering  thresh¬ 
olds  of  multiple  soil  samples.  In  the  first  set  of 
experiments  [6],  natural  peat  moss  of  negligible 
minerals  (IC  =  3.7%)  was  mixed  with  water  (m„) 
and  mineral  clay  (mcI)  to  produce  modified  soil 
samples  of  known  MC  =  ?n„/(mp  +  mct )  and 
I C  =  mcl / (mp  +  mci) .  The  modified  samples  were 
tested  in  an  insulated  box  with  the  top  open  to 
the  atmosphere,  of  a  depth  of  40  mm  and  a  cross 
section  of  90  x  90  mm  (internal  dimension).  A 
coil  heater  was  in  contact  with  the  top  surface 
for  3  min  to  initiate  the  smoldering  at  least  near 
the  coil.  In  the  second  set  of  experiments  [7],  a 
large  number  of  unmodified  natural  soil  samples 
with  different  natural  ICs  at  various  sites  of  North 
America  were  tested  to  find  MC£  and  ICc.  The 
ignition  protocol  was  different:  an  additional  layer 
of  dry  peat  of  10  mm  was  placed  between  the  coil 
heater  and  the  soil  sample. 

During  smoldering,  the  peat  is  first  dried  and 
decomposed  to  char,  and  then  char  is  oxidized 
to  ash  [2,12],  After  the  ignition  of  a  vertical  sam¬ 
ple,  a  smoldering  front  starts  to  spread  in-depth 
and  a  layer  of  ash  is  accumulated  on  the  top,  as 
illustrated  in  Fig.  1.  Depending  on  both  MC 
and  IC,  this  smoldering  front  may  become  self- 
sustaining  and  consume  most  of  the  organic  mat¬ 
ter,  or  may  not  spread  beyond  the  ignition  zone. 
As  a  first  approximation,  Frandsen’s  experiments 
can  be  modelled  as  1-D  at  the  in-depth  direction 
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flow 


z  =  CD- 


Spread 
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Emission 

A  A 

T  gases  ; 
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char  t  clay 
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wet  (peat  +  clay) 


Fig.  1.  Illustration  of  the  1-D  computational  domain  for 
in-depth  spread  of  smoldering  in  a  sample  of  peat. 


because  buoyancy  plays  a  negligible  role  in  these 
relatively  small  samples  (40  mm),  i.e.  vertical  and 
horizontal  samples  behave  the  same  way.  We  used 
the  open-source  code  Gpyro  [11]  to  implement  a 
1-D  model  and  simulate  the  ignition  and  spread 
of  smoldering  combustion  of  peat. 

2.1.  1-D  governing  equations 


The  computational  domain  is  a  40  mm  deep 
sample  with  ignition  at  the  top  free  surface  and 
insulated  at  the  bottom  (Fig.  1).  The  model  solves 
the  1-D  transient  equations  for  both  solid  and  gas 
phases  in  the  absence  of  gravity.  The  gas-phase 
temperature  is  assumed  to  be  the  same  as  the  con¬ 
densed-phase  temperature  (thermal  equilibrium). 
The  details  are  reported  in  [11];  only  the  essential 
conservation  equations  are  presented  here:  (1) 
condensed-phase  mass,  (2)  condensed-phase  spe¬ 
cies,  (3)  condensed-phase  energy,  (4)  gas-phase 
mass,  (5)  gas-phase  species,  and  (6)  gas-phase 
momentum  (Darcy’s  law).  Symbols  are  explained 
in  the  nomenclature,  and  subscripts  i,j,  and  k 
refer  to  the  number  of  condensed-phase  species, 
gas-phase  species,  and  reaction,  respectively. 


dt~  * 
djpYi) 
dt 

d{ph) d(-j dT 


dt  dz  \  dz 

d(PgJ)  dm" 
dt  +  dz 
d(Pg4'Yi)+d(m"Yj, 


fg 


d 


dt 


KdP 

m  = - 

v  dz 


dz 


dY, 


=  -xz\ 


dz 


+ <4  - 


dz 


(P  =  pgRsT) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


Each  condensed-phase  species  is  assumed  to  have 
constant  properties  (e.g.  bulk  density,  specific 
heat,  and  porosity).  All  gaseous  species  have  unit 
Schmidt  number,  and  equal  diffusion  coefficient 
and  specific  heat.  The  averaged  properties  in  each 
cell  are  calculated  by  weighting  appropriate  mass 
or  volume  fractions  [11], 

At  the  free  surface  (z  =  0),  the  convective 
boundary  condition  is  imposed  hc  o  —  10  W/m2- 
K  with  reradiation  (e  =  0.95)  and  mass  transfer 
h„,,o  =  0.02  kg/m2  ■  s  (approximation  of  Couette 
flow  [11]).  The  ambient  pressure  and  temperature 
are  assumed  to  be  atmospheric  and  300  K.  Within 
the  first  3  min,  a  heat  flux  of  30  kW/m2  is  applied 
to  simulate  the  heating  of  the  coil  heater  (the  effect 
of  the  heat  flux  level  is  explored  in  Section  4.2).  At 
the  deep  end  of  the  sample  (L  =  40  mm),  the  mass 
flux  is  set  to  zero,  and  heat  loss  is  set  with 
hc,L  =  3  W/m2  ■  K.  A  fully  implicit  formulation  is 
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adopted  for  solution  of  all  equations.  More  details 
about  the  numerical  solution  are  reported  in  [1 1]. 
Simulations  were  run  with  an  initial  cell  size  of 
Az  =  0.1  mm  (400  cells  per  domain),  and  an  initial 
time  step  of  0.01  s.  Reducing  the  cell  size  and  time 
step  by  a  factor  of  2  gives  no  significant  difference 
in  results,  so  the  calculations  are  sufficiently 
resolved. 

2.2.  Chemical  kinetics 


The  heterogeneous  reaction  in  mass  basis  is 
written  as: 

N  N 

A  +  Yy'j* gas  J  -*■  vB.kBk  +  J2vUgas  ■>’  (7) 

y=l  7=1 


where  vB,4  =  1  +  (pB/pA  -  \)xk,  and  /k  quantifies 
the  shrinkage  or  intumescence  of  the  cell  size. 
The  destruction  rate  of  condensed  species  A  in 
reaction  k  is  expressed  by  the  Arrhenius  law 


(pYA  A=)z 

Az 


Zke 


—Ek/RT 


pY A  txz 


ro,  > 


(8) 


(pYAAz)z 

(pYAAz)z  =  (pYAAz)\t=0+  f  w'JlAz(T)dT,  (9) 
Jo 


where  subscripts  d  and  /  represent  destruction  and 
formation.  The  formation  rate  of  condensed  spe¬ 
cies  B  and  all  gases  from  reaction  k  are 
and  =  (1  -  The  cor¬ 

responding  heat  of  reaction  is  Q'f  =  —co'^AHt. 

In  previous  work  [12]  the  decomposition 
schemes  with  different  complexities  were  investi¬ 
gated  using  thermogravity  (TG)  data  of  four  dif¬ 
ferent  peat  samples  from  Scotland  (SC),  Siberia 
(SI-A  and  SI-B),  and  China  (CH).  The  best  kinet¬ 
ics  scheme  was  found  to  be:  (1)  drying  (dr),  (2) 
peat  pyrolysis  (pp ),  (3)  peat  oxidation  (po ),  (4)  (3- 
char  oxidation  (/ 3o ),  and  (5)  a-char  oxidation 
(ao)  as 


Peat  •  vwdrH20  Peat  +  vWi<frH20(g) 

(dr) 

(10) 

Peat  — >  v^a-Char  +  vgtPPGas 

(PP) 

(ID 

Peat  +  vo2lB„02  ->  v^/l-Char  +  v?iPOGas  (po) 

(12) 

/1-Char  +  vo2,po02  ->  va$0  Ash  +  vgi/)oGas 

(Po) 

(13) 

a-Char  +  v02,co02  — >  vaTOAsh  +  v?i„0Gas 

(ao) 

(14) 

where  subscripts  w,p,a,f3,  and  a  represent  five 
condensed  species  (water,  peat,  a-char,  /1-char, 
and  ash),  in  addition  to  four  gaseous  species:  oxygen, 
nitrogen,  water  vapour,  and  emission  gases. 


3.  Smoldering  structure 

3.1.  Parameter  selection 

The  physical  properties  of  the  condensed- 
phase  species  are  listed  in  Table  1.  The  solid 
(ip  =  0)  physical  properties,  ps{,ksA,cPti  of  peat, 
char,  and  clay  are  selected  from  [13],  The  porosity 


Table  1 

The  physical  parameters  of  condensed-phase  species 
before  mixing  where  psi,ksi,  and  cpd  are  from  [13],  and 
pi0  is  from  [6,16], 


Species 

(i) 

Ps,i 

(kg/m3) 

Pi,  o 

(kg/m3) 

iA(,o 

(-) 

(W/m-K) 

Cp,i 

(J/kg-K) 

Water 

1000 

1000 

- 

0.6 

4186 

Peat 

1500 

110 

0.927 

1.0 

1840 

a-Char 

1300 

135 

0.896 

0.26 

1260 

/1-Char 

1300 

135 

0.896 

0.26 

1260 

Ash 

2500 

19.5 

0.992 

0.8 

880 

Clay 

2500 

1200 

0.520 

0.8 

880 

is  calculated  using  the  bulk  density  in  the  litera¬ 
ture  as  i//,.  =  (1  —  pj psi).  The  effective  thermal 
conductivity  includes  the  radiation  heat  transfer 
across  pores  as  kt  =  ksi(\  —  i/r;)  +  yoT3  where 
y=l(T4~10_3m  depends  on  the  pore  size 
(y  ~  dp  =  l/Sp)  with  a  soil  particle  surface  area 
S  ~  0.05  m2/g  [14].  The  permeability  ( K  ~  d2p) 

varies  from  10~12  to  10-9  m2  [15], 

The  sample  volume  expands  after  peat  is  mixed 
with  clay,  but  for  unsaturated  samples,  water 
occupies  the  pore  space  and  does  not  expand  the 
volume.  Thus,  the  bulk  density  of  mixed  unsatu¬ 
rated  soil  is  given  as  p  =  (1  +  MC)/[(1  —  IC)/ 
pp  +  IC /pc,\.  The  properties  of  a-char  and  /1-char 
are  assumed  to  be  the  same,  and  so  as  natural 
minerals  (ash)  and  clay.  After  combustion,  the 
residue  is  a  homogenous  mixture  of  natural  min¬ 
erals  and  clay. 

The  Scotland  (SC)  peat  sample  with  a  high 
organic  content  (IC  =  1.8%)  [12]  is  selected  as 
the  base  case.  Table  2  lists  the  kinetic  and  stoichi¬ 
ometric  parameters.  The  heat  of  oxidation  is 
related  to  the  fraction  of  oxidized  organic  matter 
and  assumed  to  be  A Hk  =  Ck(  1  —  v*).  By  integrat¬ 
ing  the  energy-flux  curve  of  TG-DSC  measure¬ 
ments  of  multiple  peat  samples  [17,18], 
Cpo  =  10  MJ/kg  and  Cm  =  Cp0  =  20  MJ/kg  are 
selected.  Under  TG  conditions  and  a  sufficient 
air  supply,  the  total  heat  of  combustion  for  a 
dry  SC  samples  is  calculated  to  be  15  MJ/kg,  sim¬ 
ilar  to  measured  values  in  [17,19],  We  relate  the 
oxygen  consumption  to  the  heat  of  oxidation 
by  assuming  a  constant  heat  of  combustion  per 
unit  of  oxygen  consumed  as  v02  k  =  A//*/(13.1 
MJ/kg)  [20], 

Table  2 


Reaction  parameters  and  gaseous  yields  of  5-step 
reactions  for  SC  peat  sample  [12], 


Parameter/k 

dr 

PP 

po 

no 

Po 

lg Ak  (lg/s-1)) 

8.12 

5.92 

6.51 

1.65 

7.04 

Ek  (kj/mol) 

67.8 

93.3 

89.8 

54.4 

112 

>h  (=1 

2.37 

1.01 

1.03 

0.54 

1.85 

vBik  (kg/kg) 

0 

0.75 

0.65 

0.03 

0.02 

A Hk  (MJ/kg) 

2.26 

0.5 

-3.54 

-19.5 

-19.5 

vo,,*  (kg/kg) 

0 

0 

0.27 

1.48 

1.49 
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Fig.  2.  Predicted  evolution  of  the  temperature  profile  at 
different  depths  for  SC  peat  with  MC  =  30%  and 
IC  =  40%. 


3.2.  Base  case 

A  successful  ignition  and  complete  combustion 
of  the  sample  is  defined  when  most  (>  95%)  of  the 
organic  matter  is  consumed,  and  only  a  small 
amount  of  char  remains  due  to  extinction  near 
the  deep  boundary.  If  ignition  succeeds,  simula¬ 
tions  show  that  the  following  spread  in  the 
40  mm  deep  samples  lasts  for  about  1  h,  similar 
to  the  observed  average  spread  rate  of  30  mm/h 
in  [6],  Computational  results  show  that  the  spread 
rate  decreases  with  MC. 

Arbitrarily,  we  choose  a  base  case  (SC  sample, 
MC  =  30%  and  IC  =  40%)  to  investigate  the  com¬ 
bustion  process  in  detail.  Figure  2  shows  the  evo¬ 
lution  of  temperature  profile.  Right  after  ignition 
(3  min),  the  smoldering  front  reaches  the  peak 
temperature  (~  900  °C),  and  then  spreads  down¬ 
wards  along  with  the  free  surface  regression.  In 
the  middle  of  the  sample  (10<z<30mm, 
t  <  25  min),  where  drying  and  peat-decomposi¬ 
tion  stages  are  important,  the  peak  temperature 
stabilizes  at  about  650  °C,  agreeing  with  the 
experimental  observation  in  [5].  Near  the  deep 
boundary,  char-oxidation  dominates  and  the  peak 
temperature  increases  upon  800  °C.  Combustion 
is  quenched  at  the  deep  end  due  to  heat  loss  and 
the  lack  of  peat. 

Figure  3a  shows  that  water  is  vaporized  within 
20  min;  peat  is  decomposed  within  27  min; 
and  char-oxidation  dominates  the  rest  of  time. 
Figure  3b  shows  that  the  peak  mass-loss  rate 
occurs  at  the  beginning  of  ignition  heating,  along 
with  a  quick  regression  of  the  fuel  bed.  Once  the 
ignition  source  is  removed,  the  mass-loss  rate  sud¬ 
denly  drops,  and  then  slowly  decreases  during  the 
peat-decomposition  stage.  The  fuel  bed  continues 
to  shrink  rapidly  due  to  both  the  mass  loss  and  the 
density  increase  from  peat  to  char.  When  the  char- 
oxidation  dominates,  the  mass-loss  rate  slightly 
decreases  due  to  both  the  end  heat  loss  and  the 
limited  oxygen  supply  with  the  accumulation  of 
mineral  residue.  The  regression  becomes  slower 


because  the  mass  loss  is  compensated  by  the  den¬ 
sity  decrease  from  char  to  ash. 

Figure  3c  shows  the  spread  at  t  =  15  min.  At 
this  instant,  the  thickness  of  smoldering  front  is 
about  25  mm,  and  there  are  three  sub-fronts  (from 
deep  to  shallow):  drying,  peat  decomposition,  and 
char  oxidation,  and  their  reaction  rates  are  on  the 
same  order  of  magnitude.  The  reaction  rate  of 
peat  oxidation  is  two  orders  of  magnitude  smaller 
than  peat  pyrolysis  (too  small  to  be  observed  in 
Fig.  3c)  because  most  of  the  oxygen  is  consumed 
in  the  char-oxidation  sub-front  and  little  oxygen 
diffuses  through  the  deeper  peat-decomposition 
sub-front.  Consequently,  very  little  of  /I-char  is 
generated,  and  the  /1-char  oxidation  is  small, 
agreeing  with  the  previous  prediction  using  a  sim¬ 
plified  plug-flow  model  [12],  Near  the  end  of  smol¬ 
dering,  only  char  is  left,  and  the  char-oxidation 
rate  is  high,  thus  increasing  the  temperatures  in 
this  region  (see  Fig.  2). 


4.  Smoldering  thresholds 

The  MC,,  measured  in  experiments  of  modified 
samples  is  plotted  against  ICc  in  Fig.  4,  and  fitted 
by  a  linear  correlation,  MCc  =  1.1  —  1.35ICc  as  in 

[6] ,  The  experimental  data  of  natural  samples  in 

[7]  is  also  plotted.  In  general,  Frandsen’s  linear 
correlation  provides  a  reasonable  reference  for 
smoldering  thresholds,  but  experimental  data 
shows  that  some  soil  samples  can  still  burn  some 
distance  above  the  line.  This  is  probably  due  to 
the  differences  in  decomposition  kinetics,  ignition 
protocol,  and  physical  properties.  The  influences 
of  these  factors  are  investigated  here  using  the  1- 
D  model. 

4. 1.  Influence  of  kinetic  parameters 

By  fixing  ICc  over  a  range  of  values,  MCc  can 
be  found  by  increasing  the  moisture  until  the 
incomplete  combustion  occurs.  On  top  of  the  SC 
peat,  the  smoldering  threshold  for  another  three 
peat  types  from  Siberia  (SI- A  and  SI-B)  and 
China  (CH)  were  also  computed  with  the  corre¬ 
sponding  kinetic  parameters  found  in  [12].  The 
SI-A  peat  also  has  a  high  organic  content 
(IC  =  2.4%),  so  the  critical  curve  similar  to  that 
of  SC  peat  is  obtained.  For  the  high-mineral  sam¬ 
ples  of  SI-B  and  CH  peat  (natural  IC  =  12. 1%  and 
18.7%),  only  one  critical  point,  MC„( natural  IC), 
is  included  in  Fig.  4. 

Figure  4  shows  that  the  computed  critical 
curve  for  both  SC  and  SI-A  samples  are  nonlin¬ 
ear.  For  SC  peat,  the  critical  curve  crosses  the 
middle  of  the  experimental  data  scatter.  The  crit¬ 
ical  curve  of  SI-A  samples  is  much  lower,  i.e. 
requires  a  much  lower  MC  to  ignite.  The  MC„ 
of  SI-B  and  CH  samples  falls  above  and  below 
the  curve  for  SC.  In  short,  the  decomposition 
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Fig.  3.  Predicted  evolution  of  (a)  mass  of  each  species  condensed  species,  (b)  total  mass-loss  rate  and  thickness  of  fuel 
bed,  and  (c)  reaction-rate  profile  at  t  =  15  min  for  SC  peat  with  MC  =  30%  and  IC  =  40%. 
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Inorganic  content,  IC  (%) 

Fig.  4.  Predicted  and  measured  critical  curves  MCc  vs.  ICc  for  the  smoldering  ignition  thresholds  of  different  soils. 
Above  the  curve,  a  soil  sample  cannot  ignite. 


Fig.  5.  Predicted  smoldering  thresholds  with  different  ignition  protocols. 


2680 


X.  Huang  et  al.  I  Proceedings  of  the  Combustion  Institute  35  (2015)  2673-2681 


Inorganic  content,  IC  (%) 


Fig.  6.  Predicted  smoldering  thresholds  for  a  combination  of  pp,ksp,  and  Cm,  varying  in  a  wide  range. 


kinetics  have  a  significant  influence  on  smoldering 
thresholds. 

4.2.  Influence  of  ignition  protocol 

The  influence  of  ignition  protocol  on  the  smol¬ 
dering  thresholds  is  investigated  by  varying  the 
external  heat  flux  in  the  range  from  15  to 
30  kW/m2.  The  results  are  shown  in  Fig.  5.  The 
minimum  heat  flux  to  ignite  a  dry  SC  peat,  whose 
MC  is  in  equilibrium  with  the  ambient  (MC,/r 
~  10%  [17,18]),  is  found  to  be  10.5  kW/m2  for 
3  min  of  heat,  and  5.8  kW/m2  for  30  min.  This 
critical  heat  flux  is  much  lower  than  the  typical 
value  for  flaming  fires  (~30kW/m2)  [21].  The 
dry-peat  ignition  protocol  used  in  the  second  set 
of  experiments  of  [7]  is  also  simulated  by  adding 
1  cm  of  dry  peat  on  the  top  which  is  heated  for 
the  same  3  min  under  30  kW/m2. 

For  organic  soils  with  IC  <  30%,  increasing  the 
heat  flux  extends  the  smoldering  threshold 
because  the  drying  of  the  whole  sample  bed 
induced  by  the  ignition  source  increases.  Slightly 
above  the  critical  curve  (e.g.  point  A  — >  A1),  the 
reaction  zone  is  quickly  quenched  (<5  min)  near 
the  top  surface  within  the  preheated  zone.  There¬ 
fore,  for  these  cases,  the  threshold  is  defined  as  the 
critical  moisture  of  ignition  (MCCiig).  For  mineral 
soils  with  IC  >  30%,  the  influence  of  the  ignition 
protocol  becomes  negligible.  Slightly  beyond  the 
critical  curve  (e.g.  B  — >  B'),  the  smoldering  front 
is  able  to  spread  out  of  the  preheated  zone  for  a 
distance  before  extinction.  These  cases  were 
recorded  as  “partial  burn”  in  the  experiments  of 


[6,7],  therefore,  the  threshold  is  defined  as  the  crit¬ 
ical  moisture  of  extinction  (MCtJ. 

The  results  show  that  the  dry-peat  ignition 
protocol  is  nearly  equivalent  to  a  heat  flux  of 
55  kW/m2  without  the  additional  layer  of  dry 
peat.  Once  ignited,  a  sustained  smoldering  front 
is  generated  so  that  the  threshold  is  for  extinction. 
When  the  heat  flux  is  very  low  (~15  kW/m2),  the 
whole  critical  curve  declines,  and  the  threshold  is 
for  ignition. 

4.3.  Sensitivity  to  properties  and  heat  of  combustion 

The  physical  properties  of  peat  are  expected  to 
vary  somehow  for  different  ecosystems.  This 
might  affect  the  smoldering  thresholds.  Here,  the 
sensitivity  of  the  thresholds  is  investigated  with 
combinations  of  the  three  most  important  param¬ 
eters  varying  over  wide  ranges  of  values  found  in 
the  literature:  bulk  density  of  peat  (pp  =  90- 
130  kg/m3  [6]),  solid  thermal  conductivity  of  peat 
(ksp  =  0.8—1 .2  W/m  K  [13]),  and  heat  of  char  oxi¬ 
dation  (Cao  =  Cfj0  =  15-25  MJ/kg  [17]).  Figure  6 
reveals  that  the  range  of  predicted  smoldering 
thresholds  covers  most  of  the  experimental  data. 

Computational  results  show  that  MCc 
increases  monotonously  with  the  heat  of  combus¬ 
tion,  but  decreases  with  the  bulk  density  and  the 
thermal  conductivity.  Therefore,  the  case  with 
maximum  Cao  and  minimum  pp  and  ksp  gives 
the  upper  boundary  (see  Fig.  6).  As  Cao  increases, 
the  heat-generation  rate  increases,  overcoming  the 
heat-sink  effect  of  water  and  minerals.  The 
increase  in  ksp  is  almost  equivalent  to  the  increase 
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in  IC  regarding  the  threshold,  thus  lowering  the 
critical  curve  (i.e.  more  difficult  to  smolder),  espe¬ 
cially  at  low  IC  values.  As  pp  increases,  the  mass 
concentrations  of  water  and  clay  also  increase 
under  the  same  MC  and  IC.  Therefore,  the 
increasing  heat-sink  effect  moves  the  critical  curve 
towards  lower  MC  values,  agreeing  with  experi¬ 
mental  observations  in  [22]. 

5.  Conclusions 

In  this  work,  a  comprehensive  1-D  model  of  a 
reactive  porous  media  is  implemented  in  the  open- 
source  code  Gpyro  to  investigate  the  smoldering 
combustion  of  peat  with  a  5-step  (including  dry¬ 
ing)  heterogeneous  kinetics.  Two  sets  of  small- 
scale  experiments  [7,6]  are  simulated  for  the  first 
time,  and  the  transient  temperature,  species,  reac¬ 
tion  profiles,  and  surface  regression  of  the  fire  are 
studied.  The  predicted  smoldering  thresholds 
related  to  the  critical  moisture  and  inorganic  con¬ 
tents  are  nonlinear,  as  opposed  to  previously 
reported  linear  correlation,  and  show  a  better 
agreement  with  the  experimental  results  for  a  wide 
range  of  soil  types.  The  smoldering  thresholds  are 
found  to  depend  on  the  decomposition  kinetics, 
physical  properties,  and  the  ignition  protocol. 
The  results  reveal  that  the  threshold  values  found 
by  Frandsen  for  organic  soils  are  due  to  ignition, 
and  for  mineral  soils  are  due  to  extinction.  This  is 
the  first  time  that  a  physics-based  model  of  smol¬ 
dering  peat  fires  is  developed,  thus  helping  to 
understand  this  important  natural  and  widespread 
phenomenon. 
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